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Using events in which one of two neutral B mesons from the decay of an T(45) meson is fully 
reconstructed, we determine parameters governing decay (ATd/Td), CP and T violation {\q/p\), 
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and CP and CPT violation (Rez, Imz). The results, obtained from an analysis of 88 million T(4S') 
decays recorded by BABAR,, are 

sgn(ReAcp)Ar<i/rd = -0.008 ±0.037 (stat.) ±0.018 (syst.) [ -0.084 , 0.068 ] , 

\q/v\ = 1.029 ±0.013 (stat.) ±0.011 (syst.) [ 1.001,1.057], 

(ReAcp/|Acp|)Rez = 0.014 ±0.035 (stat.) ± 0.034 (syst.) [-0.072,0.101], 

Imz = 0.038 ±0.029 (stat.) ±0.025 (syst.) [ -0.028 , 0.104 ] . 

The values inside square brackets indicate the 90% confidence-level intervals. These results are 
consistent with Standard Model expectations. 

PACS numbers: 13.25.Hw, 12.15.Ff, ll.30.Er 



In this Letter, we provide a direct limit on the total 
decay-rate difference AF^ between the Bd mass eigen- 
states, and set limits on CP, T, and CPT violation inher- 
ent in the mixing of neutral B mesons. In the Standard 
Model CPT violation is forbidden and the other effects 
are expected to be non-zero but smalL but new physics 
could provide enhancements P, Q S 13 ■ We test these 
predictions by analyzing the time dependence of decays of 
the T{AS) resonance in which one neutral B meson (Brcc) 
is fully reconstructed and the flavor of the other B (Stag) 
is identified as being either or 5° . The Brec sample is 
composed of flavor- and CP-eigenstate subsamples, iJfiav 
and Bcp- We reconstruct the flavor eigenstates -Bflav 
= D(*)-7r+(p+,a;i^) and J/iIjK*°{^ K+tt-) and the CP 
eigenstates Bcp =J/^pK°, iP{2S)K°, XciK°, and J/^pK1. 
The flavor of the B that is not completely reconstructed 
is "tagged" on the basis of the charges of leptons and 



kaons, and other indicators 



The data come from 



million TUS) BB decays collected with the BABAR 
detector [7| at the PEP-II asymmetric-energy B Factory 
at SLAG. 

The light and heavy Bd mass eigenstates Bl h are su- 
perpositions of B° and B°. This mixing is a consequence 
of transitions between B" and B" through intermediate 
states. Flavor oscillations between B° and B° occur with 
a frequency Amd = niH — 'mL- A state that is initially B° 
(B") will develop a B" (B") component over time, whose 
amplitude is proportional to a complex factor denoted 
q/p (p/q) [ll- Since \q/p\ ~ 1 in the Standard Model, 
this factor is usually assumed to be a pure phase. 

The most general time dependence allowed for the de- 
cays of the two neutral B mesons coming from an T{AS) 



ays ( 



diV 
dAt 

2 



oc e 



|a_r . f ATdAt 



■ cosh 



+ 



\a+ 



■ cos(AmdAi) — Re{a*^_a^) sinh 



/ AFrfAt 



Im(a!i_a_) sin(ATOc;Ai) 



(1) 



where At = tree — ttag is the signed difference in proper 
decay times, F^ is the mean decay rate of the two neutral 
mass eigenstates, and AF^ = Th — F^, is their decay- 



rate difference. The values of the complex parameters 
a± differ for the various combinations of flavor and CP 
eigenstates into which the B mesons decay 6] . 

In the simplest picture, where AF^ = 0, and CP, T, 
and CPT violation in mixing are neglected, if the fully 
reconstructed state is a flavor eigenstate the time distri- 
butions dN/dAt with perfect tagging are proportional 
to e^'"''!'^*! (1 ± cos{AmdAt)) . In practice, the tagging is 
imperfect and its performance is measured directly from 
the data. Imperfect tagging reduces the coefficient of 
cos{AmdAt) by a factor 1 — 2u' called the dilution, where 
w is the probability of tagging incorrectly. 

B decays to a CP eigenstate fcp are conveniently pa- 
rameterized by Xcp = (q/p) (Acp/Acp), where Acp 
(Acp) is the amplitude for B° fcp {B° fcp). CP 
violation is characterized by Xcp 7^ Tjcp where r]cp = ±1 
is the final state's CP eigenvalue. The CP violation ob- 
served in decays like B J/ip K° 0] involves inter- 
ference between decays with and without net oscillation, 
and leads to ImAcp 7^ 0. Other possible sources of CP 
violation are \q/p\ ^ 1 and |^cp/-4cp| 7^ 1. We include 
a test of the former possibility here. 

The time distributions dN/dAt for the Bcp samples, 
in the simplest picture (defined above) and with perfect 
tagging, are proportional to 

g-r,|At| ^ |^^^|2 ^ _ \Xcp\^)cos{AmdAt) 

T 2ImAcpsin(AmdAt)] . (2) 

In the Standard Model we have Xcp — — e^^*^ for 
J/i/jK'^ with the approximation AF^ = 0, where /3 = 
arg f— y,^V*K/Vt^V*K] is one of the angles of the trian- 
gle [ll| that represents the unitarity of the quark mixing 
matrix Vy. Since \Xcp\ = 1, the cos(AmrfAt) term is ab- 
sent. Again, wrongly tagged events reduce the amplitude 
of the oscillatory terms. 

To measure AF^, or CP, T, or CPT violation in mix- 
ing alone we need to find small deviations from these 
simple patterns. Other effects that can mimic the be- 
havior we seek must be included in the analysis. Among 
these are asymmetries in the response of the detector to 
B" and B*^ decays and interference between dominant 
and suppressed decay amplitudes to flavor eigenstates, 
both those that are fully reconstructed and those that 
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contribute to tagging 0| . 

The time dependence of the Bcp sample includes a 
sinh(Ar(iAt/2) term that is effectively linear in AF^;, 
while the flavor sample has an effective second-order sen- 
sitivity to ATd through a cosh(ArdAt/2) term. Un- 
tagged data are included in this analysis and improve 
our sensitivity to AF^ since the contributions of AF^- 
dependent terms do not depend on whether Stag is a 
_B" or i?'^. With our sample sizes and small measured 
value of ATd, the Bcp sample dominates our deter- 
mination of AFrf/Frf. While Anid has been well mea- 
sured previously there is only a weak hmit, 
|AFrf|/Fd < 0.18 at 95% CL [l^, on AF^. A recent the- 
oretical calculation gives AFd/r^ = —0.003 0. 

Violation of CP and T in mixing leads to a difference 
between the — > B° and i?" B^ transition rates 
proportional to — 1- Our sensitivity to \q/p\ comes 

mostly from the large flavor-eigenstate sample. Previ- 
ous measurements, obtained assuming AF^ = 0, give 
\q/p\-l = (-0.7±6.4)xl0-3 Il7|. The Standard Model 
expectation is \q/p\ - 1 ^ (2.5 6.5) x 10"^ ,4 

CPT violation in mixing enters the time dependence 
through the complex quantity 

Arud - ^ATd ' ^ ' 

where 6md {6Td) is the B^ — B^ difference of effective 
mass (decay rate) expectation values for the B^ and B^ 
flavor eigenstates. A non-zero value of either 5md or 
5Td is only possible if both CP and CPT are violated. 
The dominant contribution of Imz to the time depen- 
dence is through the coefficient of sm.{AmdAt) for fla- 
vor eigenstates, while Rez contributes primarily to the 
coefficients of cosh(AFdAt/2) « 1 and cos{AmdAt) for 
CP eigenstates. The measurement of z presented here is 
more general than previous analyses based on B decays, 
which obtained Imz = 0.040 ± 0.032 ± 0.012 [3, and 
Rez = 0.00±0.12±0.02, Imz = -0.03 ± O.Oljt 0.03 [13, 
and complements earlier limits on the — K'^ mass dif- 
ference SmK/rriK < 10~^^ @. 

Interference effects between the amplitudes for domi- 
nant decays of flavor eigenstates (e.g., B° — > D~tt~^) and 
for doubly-CKM-suppressed (DCS) decays (e.g., B^ 
D~'k'^) are analogous to the interference familiar in de- 
cays to CP eigenstates 0. They thus affect, in partic- 
ular, the sin(ATOc(Ai) terms and have the potential to 
obscure a similar contribution from Imz. The size of the 
DCS interference relative to the dominant B^ decay is 
governed by As/ and Xst, for Bflav and i?tag states, re- 
spectively. These parameters are defined analogously to 
Acp, and we expect \\Bf,Bt\ ~ k/p|lKbKyK*bK.dl - 
0.02\q/p\ g. There are similar interference contribu- 
tions from DCS amplitudes for decays, governed by 
A^^ and Xg^. We write X-Bf^m = ^/^'BfM so l^Bf^ml ~ 
0.02|p/q|. The iJflav and Stag samples are ensembles of fi- 
nal states that each contribute to the expected decay-rate 



distributions with different amplitudes. We find that, 
working to first order in the small quantities \XBf,Bt\, 
l^BfBtl^ 1^1' ^^'^ \l/p\ ^ 1' the cumulative effect of each 
ensemble does not modify the expected decay-rate dis- 
tributions, once XBf.Bt and X-^j^^ are reinterpreted as 
effective parameters. 

We combine all the data for the CP eigenstates, taking 
into account the CP eigenvalue of the final state. We as- 
sume |^cp/-4cp| = 1 (but vary this ratio as a systematic 
study) as expected theoretically at the 10~^ level jT^ and 
as supported by the average of B-Factory measurements 
of states of charmonium and or , for which it was 
found \Acp/Acp\ = 0.949 ±0.045 [i[ig, when AF<j = 
and \q/p\ = 1 are assumed. 

The time interval At between the two B decays is cal- 
culated from the measured separation Az between the de- 
cay vertices of i?rec and Btag along the collision axis 0| . 
We determine the position of the i?rec vertex from its 
charged tracks. The Btag decay vertex is determined 
by fitting to a common vertex tracks not belonging to 
the i3rcc candidate, employing constraints from the beam 
spot location and the B^ec momentum 0. The r.m.s. At 
resolution for 99.7% of the events used is 1.0 ps, to be 
compared with {\At\) ~ 1.5 ps. 

Resolution effects for signal events are modeled by con- 
volving the idealized decay rate with a sum of three Gaus- 
sian distributions, two of whose widths and biases are 
scaled with each event's estimated At uncertainty (ta*- 

We use four mutually-exclusive categories to assign 
tags, based on kinematic, particle type, and charge infor- 
mation . There are separate reconstruction efficiencies 
and mistag probabilities for B^ and B^ tags, to accom- 
modate differences in detector response to B^ and B^ 
decays. In addition, we introduce a linear dependence of 
the mistag probability on cta*, except for events tagged 
with a high-momentum lepton. 

Backgrounds are primarily due to misreconstructed 
Broc candidates and are studied in data using mass or 
energy sidebands. Events are assigned signal and back- 
ground probabilities based on their proximity to the sig- 
nal peak. We model backgrounds with empirical At dis- 
tributions that can accommodate contributions from de- 
cays with a range of lifetimes. 

The parameters of primary interest in this analysis 
are sgn(Re Acp)AFd/Fd, \q/p\, (Re Acp/|Acp|) Re z, and 
Imz. We cannot determine AVd/Td and Rez directly be- 
cause both occur multiplied by Re Acp in their dominant 
contributions to the decay rate. They are thus subject 
to a sign ambiguity, which can be resolved by relying on 
additional information from the unitarity trian^. The 
average lifetime t^o = 1/Td is fixed at 1.542 ps The 
parameters Arud and ImAcp/jAcpl are determined to- 
gether with the main parameters as cross checks against 
earlier measurements . The terms proportional to 

the real parts of the effective DCS parameters XBf,Bt and 
Xgr-g^, are small and therefore neglected in the nomi- 
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TABLE I: Fit results allowing (z free) or not allowing (z = 0) 
CPT violation in B°B° oscillations. 



Parameter 


z Free 


z = 


sgn(ReAcp)Ard/rd 


-0.008 ±0.037 


-0.009 ± 0.037 




1.029 ±0.013 


1.029 ±0.013 


(ReAcp/|Acp|)Rez 


0.014 ± 0.035 




Imz 


0.038 ± 0.029 





nal fit model, while the imaginary parts and magnitudes 
of these effective parameters are treated as independent 
variables. For all sets of nonleptonic flavor eigenstates 
analyzed, the magnitude of each |A| is fixed to 0.02 (up to 
a factor \q/p\ or \p/q\) but Im A/|A| is left unconstrained. 
The decay model uses 26 more parameters to model the 
effects of experimental At resolution (10), B'^ /B'^ tagging 
capabihty (11), and reconstruction and tagging efficien- 
cies (5). An additional 22 parameters model the levels 
and Ai dependence of backgrounds. A total of 58 free pa- 
rameters are determined with a simultaneous unbinned 
maximum-likelihood fit to the At distributions of CP and 
flavor-eigenstate samples . 

Table n summarizes the results of fits allowing (z free) 
or not allowing (z = 0) CPT violation in B^B^ oscilla- 
tions. The largest statistical correlations involving the 
parameters of interest are between \q/p\ and parameters 
modeling B'^B^ asymmetries in reconstruction efficiency 
and mistag probabilities, and between Imz and the DCS 
contributions to i?tag decay amplitudes. The fitted val- 
ues of Am^ and Im Acp / 1 Aopjare consistent with recent 
_B-Factory measurements Isl, [lOl, [ij, [iSj ■ When z is fixed, 
the value of ImA(7p/|Acp| decreases by 0.011, equal to 
15% of the statistical uncertainty on Im A(7p/|Acp| which 
is consistent with the correlations observed in the fit with 
z free, while the value of and uncertainty in Am^ are un- 
changed. No statistically significant B^-B'^ differences in 
reconstruction and tagging efficiencies are observed. 

We have used data and Monte Carlo samples to vali- 
date our analysis technique. Tests with large, parameter- 
ized Monte Carlo samples demonstrate that the observed 
statistical uncertainties and correlations are consistent 
with expectations. Analyses of Monte Carlo samples gen- 
erated with a detailed detector simulation verify that the 
analysis procedure is unbiased. Fits to data subsamples 
selected by tagging category, running period, and B^^c 
decay mode give consistent results. Changes to the algo- 
rithms used to estimate At and (jAt or to their allowed 
ranges also have no statistically significant effect. Fits to 
samples of charged B decays, in which no oscillations are 
present, give the expected results. 

We identify four general sources of systematic uncer- 
tainty with the contributions shown in Table |n] for the 
fit in which z is free [^. The first is possible bias in 
the event selection and fit method: we see no evidence 
of such bias when analyzing Monte Carlo samples and 



TABLE II; Summary of systematic uncertainties (z free). 





sgn(Re Acp) 




Re Acp/ Acp 




Source 


xAFd/Fd 


k/pl 


X Rez 


Imz 


Analysis Method 


0.006 


0.007 


0.005 


0.016 


At Resolution 


0.013 


0.003 


0.008 


0.016 


Signal Properties 


0.010 


0.008 


0.033 


0.009 


BG Properties 


0.005 


0.003 


0.007 


0.004 


Total 


0.018 


0.011 


0.034 


0.025 



assign the statistical uncertainty of these checks as a sys- 
tematic uncertainty on the final results. The second is the 
At measurement. The choice of parameterization of the 
resolution function dominates this uncertainty, but as- 
sumptions about the beam spot and detector alignment 
contribute as well. Assumptions about the properties of 
signal T(4S') Srcc-Btag decays include the values of the 
lifetime, \Acp / Acp\^ and DCS parameters, and are the 
third source of systematic uncertainty. Uncertainties in 
the size and At distributions of background events (BG) 
incorrectly identified as T{AS) Broc-Btag make small 
contributions to the systematic uncertainties. 

Different sources dominate the systematic uncertainty 
for each parameter. Most systematic uncertainties are 
determined with data and will decrease with addi- 
tional statistics. The largest single source of uncer- 
tainty is the contribution of the DCS parameters to 
(Re Acp /I Acp I) Rez and it is estimated by varying the 
DCS phase parameters over their full allowed range, and 
\ABf /Asfl and \A-gj/A-gj:\ over the range 0-0.04. Sys- 
tematic uncertainties on sgn(Re Acp)Arc(/rrf and \q/p\ 
for the analysis assuming z = were evaluated similarly 
as ±0.018 and ±0.011, respectively. 

Using the world-average value of Am^ [^, we derive 
the value sgn(Re Acp)Ard/ATO<i = -0.011 ±0.049(stat.) 
±0.024(syst.), corresponding to the range [-0.112,0.091] 
at the 90% confidence level, from the fit results with 
z free. The limit on CP and T violation in oscilla- 
tions is independent of and consistent with our previous 
measurement based on an analysis of inclusive dilepton 
events 2(1. Using Eq. © and taking the world-average 
Bd mass Q, we derive \6md\/mBa < 1-0 x 10"^'' and 
-0.156 < STd/Td < 0.042 at the 90% confidence level. 
Figure Q shows the results of the fit with z free in the 
{\q/p\ — 1, |z|) plane, compared to the previous BABAR 
measurement of \q/p\, and to Standard Model expecta- 
tions. 

Conventional analyses of oscillations and CP violation 
in the Bd system neglect possible contributions from sev- 
eral sources that are expected to be small in the Stan- 
dard Model. This analysis includes these effects and 
finds results consistent with Standard Model expecta- 
tions. While the Standard Model predictions for AF^, 
\q/p\, and z are well below our current sensitivity, higher- 
precision measurements may still bring surprises. 
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FIG. 1: Favored regions at 68% confidence level in the {\q/p\ — 
1, |z|) plane determined by this analysis and by the BABAR 
measurement of the dilepton asymmetry |20|| . Labels reflect 
the requirements that both CP and T be violated if \q/p\ 7^ 
1 and that both CP and CPT be violated if |z| / 0. The 
dilepton measurement constrains \q/p\ without assumptions 
on the value of |z|. The Standard Model expectation of \q/p\ — 
1 = (2.5 - 6.5) X 10-"* is obtained from Ref. 0. 
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